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A search for fermionic top quark partners T of charge 2/3 is presented. The search is 
carried out in proton-proton collisions corresponding to an integrated luminosity of 
19.7fb^^ collected at a center-of-mass energy of = 8 TeV with the CMS detector at 
the LHC. The T quarks are assumed to be produced strongly in pairs and can decay 
into tH, tZ, and bW. The search is performed in five exclusive channels: a single¬ 
lepton channel, a multilepton channel, two all-hadronic channels optimized either 
for the bW or the tH decay, and one channel in which the Higgs boson decays into 
two photons. The results are found to be compatible with the standard model expec¬ 
tations in all the investigated final states. A statistical combination of these results is 
performed and lower limits on the T quark mass are set. Depending on the branch¬ 
ing fractions, lower mass limits between 720 and 920 GeV at 95% confidence level are 
found. These are among the strongest limits on vector-like T quarks obtained to date. 
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1 Introduction 

The discovery of a Higgs boson with a mass of 125 GeV by the ATLAS [1J and CMS lU |3 col¬ 
laborations motivates the search for exotic states involving the newly discovered particle. The 
nature of electroweak symmetry breaking and the mechanism that stabilizes the mass of the 
Higgs particle are not entirely clear. These questions could be explained by physics beyond the 
standard model (SM), such as supersymmetry. Non-supersymmetric explanations are given 
by little Higgs models ||11|5|, models with extra dimensions |0|Zl, and composite Higgs mod¬ 
els ISl-El in which the Higgs boson appears as a pseudo-Nambu-Goldstone boson [j9l. These 
theories predict the existence of heavy vector-like quarks. The left-handed and right-handed 
components of vector-like quarks transform in the same way under the electroweak symmetry 
group, in contrast to the SM fermions, which transform as chiral particles under the SM sym¬ 
metry group SU(3)c X SU(2)l x U(1)y. This property of the vector-like quarks allows direct 
mass terms in the Lagrangian of the form mxpxp that do not violate gauge invariance. As a con¬ 
sequence, and in contrast to the other quark families, vector-like quarks do not acquire their 
mass via Yukawa couplings. In many of the models mentioned above the vector-like quarks 
couple predominantly to the third generation quarks only. This means that they may have 
the following three decay modes: tH, tZ, and bW KlOll . A model of vector-like T quarks with 
charge 2/3 e, which are produced in pairs via strong interaction, is used as a benchmark for this 
analysis. 

A fourth generation of chiral fermions, replicating one of the three generations of the SM with 
identical quantum numbers, is disfavored by electroweak fits within the framework of the 
SM mn . This is mostly because of large modifications of the Higgs production cross sections 
and branching fractions {B), if a single SM-like Higgs doublet is assumed. Heavy vector-like 
quarks decouple from low energy loop-level electroweak corrections and are not similarly con¬ 
strained by the measurements of the Higgs boson properties UTOl . 

Early T quark searches by the GMS Gollaboration ||12HT4| have assumed 100% branching frac¬ 
tions to various final states. More recent searches HTSl do not make specific assumptions for the 
branching fractions. Searches for T quarks have been performed also by the ATLAS Gollabo¬ 
ration, setting lower limits on the T quark mass ranging from 715 to 950 GeV, for different T 
quark branching fractions IfTOtiTSlI . 

In this paper, results of searches for T quark production in proton-proton collisions, using the 
GMS detector at the GERN LHG, are presented for five different decay modes. One of the 
searches is inclusive and sets limits for all possible branching fractions. This analysis is 
based on leptonic final states and is described in Section |5.1| The other four analyses have a 
good sensitivity in optimized regions, but they do not cover the full range of branching frac¬ 
tions. The analysis described in Section |5.2| is specifically optimized to find T —^ bW decays. 
The searches presented in Section [53| and Section [5A| are optimized for all-hadronic final states 
in the decays T —> bW and T —)■ tH. The search discussed in Section [53| is sensitive to T —> tH 
decays, where the Higgs boson decays to a pair of photons. The two analyses presented in Sec¬ 
tions |5^ and |5^ are discussed in detail in separate publications llTSlIl^ . The remaining three 
analysis are published here for the first time. 

The GMS detector is briefly described in Section Section describes the data and the simu¬ 
lated samples. Section|^gives details about the reconstruction techniques used by the analyses. 
Section describes the combination and the treatment of systematic uncertainties. Section 
presents the results of the combination. 
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3 Event samples 


2 The CMS detector 

The central feature of the CMS apparatus is a superconducting solenoid of 6 m internal diame¬ 
ter, providing a magnetic field of 3.8 T. Within the solenoid volume are a silicon pixel and strip 
tracker, a lead tungstate crystal electromagnetic calorimeter (ECAL), and a brass and scintilla¬ 
tor hadron calorimeter (HCAL), each composed of a barrel and two endcap sections. Muons 
are measured in gas-ionization detectors embedded in the steel flux-return yoke outside the 
solenoid. Extensive forward calorimetry complements the coverage provided by the barrel 
and endcap detectors. 

In the region of pseudorapidity \rj\ < 1.74 Il20ll . the HCAL cells have widths of 0.087 in r] and 
0.087 radians in azimuth {(p). In the f]-(p plane, and for \f]\ < 1.48, the HCAL cells map on to 
5x5 ECAL crystals arrays to form calorimeter towers projecting radially outwards from close 
to the nominal interaction point. At larger values of \tj\, the size of the towers increases and the 
matching ECAL arrays contain fewer crystals. Within each tower, the energy deposits in ECAL 
and HCAL cells are summed to define the calorimeter tower energies, subsequently used to 
provide the energies and directions of hadronic jets. 

The electron momentum is estimated by combining the energy measurement in the ECAL with 
the momentum measurement in the tracker. The momentum resolution for electrons with 
transverse momentum pj 45GeV from Z —> ee decays ranges from 1.7% for nonshowering 
electrons in the barrel region to 4.5% for showering electrons in the endcaps ||2T1. The energy 
resolution for photons with transverse energy Ej ^ 60GeV varies between 1.1% and 2.6% in 
the EGAL barrel, and from 2.2% to 5% in the endcaps Il22l . 

The silicon tracker measures charged particles within the pseudorapidity range \p\ < 2.5. It 
consists of 1440 silicon pixel and 15148 silicon strip detector modules. Lor nonisolated particles 
of 1 < Pt < lOGeV and \r]\ < 1.4, the track resolutions are typically 1.5% in pj and 25-90 (45- 
150) pm in the transverse (longitudinal) impact parameter II23II . 

A more detailed description of the GMS detector, together with a definition of the coordinate 
system used and the relevant kinematic variables, can be found in Ref. EOl . 


3 Event samples 


This analysis makes use of data recorded with the GMS detector in proton-proton collisions at 
a center-of-mass energy of v^s = 8 TeV corresponding to an integrated luminosity of 19.5 fb~^ 
for the analysis described in Section 5.1 and 19.7fb^^ for the other analyses. 


Events are selected by a multi-stage trigger system. The single-lepton channels are based on 
single-muon and single-electron triggers. The single-muon sample is obtained by the require¬ 
ment of an isolated muon candidate, with high-level trigger thresholds of pj > 24 GeV (in¬ 
clusive search. Section [5T] | or pj > 40 GeV (single-lepton search. Section [5^ . In the electron 
sample, a single isolated electron trigger with pj > 27 GeV is required. Multilepton events are 
selected by requiring at least two lepton candidates, one with pj > 17 GeV and the other with 
Pt > 8 GeV in the high-level trigger. The all-hadronic final states require large hadronic activ¬ 
ity in the detector, namely that the scalar pj sum of reconstructed jets is larger than 750 GeV. 
This quantity is evaluated in the high-level trigger from jets with pj > 40 GeV using calorime¬ 
ter information only. Lor searches in the diphoton final state, two photons are required. The 
photon Ej thresholds in the high-level trigger are 26 (18) GeV and 36 (22) GeV on the leading 
(subleading) photon, depending on the running period. 
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The contributions from SM processes are generally predicted using simulated event samples. 
For some backgrounds, however, the simulations are not fully reliable, and control samples of 
data are used to determine their contribution. The background estimation for the individual 
channels is discussed in Section |3 

Standard Model background events are simulated using powheg vl.O ||24]-(26l for tf and single t 
production; MadGraph 5.1 127] for W+jets, Z+jets, tfW, and ttZ production; and PYTHIA 6.426 
[2?| for WW, WZ, ZZ, and tfH processes. 

For W+jets and Z+jets production, samples with up to four partons are generated and merged 
using the MLM scheme with kj jets 112^13011 . The CTEQ6M parton distribution functions (PDF) 
are used for POWHEG, while for the other generators the CTEQ6L1 11311 PDEs are used. In all 
cases, PYTHIA 6.426 Il28l is used to simulate the hadronization and the parton showering. 

The TT signal process is simulated using MadGraph 5.1, allowing up to two additional hard 
partons. A series of mass hypotheses between 500 and 1000 GeV are generated in steps of 
100 GeV. The inclusive cross sections for the signal samples and the tt samples are calculated 
at next-to-next-to-leading order (NNLO) for gg —> tf -|- X. The fixed-order calculations are 
supplemented with soft-gluon resummations having next-to-next-to-leading logarithmic accu¬ 
racy Il32l . The tf cross sections are computed based on the TOP++ v2.0 implementation using 
the MSTW2008nnlo68cl PDFs and the 5.9.0 version of LHAPDF Il32ll33l . The tf cross section 
is computed to be 252.9 pb, assuming a top quark mass of 172.5 GeV. The model-independent 
cross sections calculated for the signal samples are listed in Table 

Table 1: The NNLO TT pair production cross section for different values of the T quark mass. 


T quark mass 
(GeV) 

Production 
cross section (pb) 

500 

0.59 

600 

0.17 

700 

0.059 

800 

0.021 

900 

0.0083 

1000 

0.0034 


Minimum bias interactions are generated using PYTHIA and are superimposed on the simu¬ 
lated events to mimic the effect of additional proton-proton collisions within a single bunch 
crossing (pileup). The pileup distributions of the simulated signal and background events 
match that observed in data, with an average of 21 reconstructed collisions per beam crossing. 


4 Event reconstruction 

Tracks are reconstructed using an iterative tracking procedure 11231 . The primary vertices are 
reconstructed with a deterministic annealing method 1341 from all tracks in the event that are 
compatible with the location of the proton-proton interaction region. The vertex with the high¬ 
est is defined as the primary interaction vertex (PV), whose position is determined 

from an adaptive vertex fit 135 1. 

The particle-flow event reconstruction algorithm l36l l37l reconstructs and identifies each indi¬ 
vidual particle, using an optimized combination of information from the various elements of 
the GMS detector. The energy of muons is obtained from the curvature of the corresponding 
track. The energy of electrons is determined from a combination of the electron momentum at 
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4 Event reconstruction 


the PV as determined by the tracker, the energy of the corresponding ECAL cluster, and the en¬ 
ergy sum of all bremsstrahlung photons spatially compatible with originating from the electron 
track. The energy of charged hadrons is determined from a combination of their momentum 
measured in the tracker and the matching ECAL and HCAL energy deposits, corrected for zero 
suppression effects and for the response function of the calorimeters to hadronic showers. Ei- 
nally, the energy of neutral hadrons is obtained from the corresponding corrected ECAL and 
HCAL energies. 

Muon (electron) candidates are required to originate from the PV and to be isolated within 
AR = + (A<^)2 < 0.4 (0.3) around the lepton direction, where A?/ (A^) indicates the 

difference in pseudorapidity r] {(p) from the lepton direction. The degree of isolation is quanti¬ 
fied by the ratio of the pj sum of all additional particles reconstructed in the isolation cone to 
the Pt of the lepton candidate. This ratio for a muon (electron) is required to be less than 0.12 
(0.10). Together with the lepton identification requirements, the isolation conditions strongly 
suppress backgrounds from jets containing leptons. 

Photons are identified as ECAL energy clusters not linked to the extrapolation of any charged 
particle trajectory to the ECAL. The energy of photons is directly obtained from the ECAL 
measurement, corrected for zero-suppression effects. In the ECAL barrel section, an energy 
resolution of about 1% is achieved for unconverted or late-converting photons in the tens of 
GeV energy range. The remaining barrel photons are measured with an energy resolution of 
about 1.3% up to |//| = 1, rising to about 2.5% at |/y| = 1.4. In the endcaps, the resolution of 
unconverted or late-converting photons is about 2.5%, while all other photons have a resolution 
between 3 and 4% Il38] . 

Eor each event, hadronic jets are reconstructed by applying the anti-kj (AK) algorithm Il3^l40l 
and/or the Cambridge-Aachen (CA) i4TI jet clustering algorithms to the reconstructed parti¬ 
cles. The AK algorithm is used with a jet size parameter of 0.5 (AK5 jets). In some analyses 
both algorithms are used. The algorithms are applied independently of each other to the full 
set of reconstructed particles. Charged particles that do not originate from the PV are removed 
from the jets. The momentum of each jet is defined as the vector sum of all particle momenta 
in the jet cluster, and is found in the simulation to be within 5% to 10% of the true particle-level 
momentum over the whole pj spectrum and detector acceptance. Jet energy corrections are de¬ 
rived from the simulation, and are confirmed with measurements of the energy balance of dijet 
and photon+jet events Il42ll . The jet energy resolution is typically 15% at 10 GeV, 8% at 100 GeV, 
and 4% at 1 TeV, to be compared to about 40%, 12%, and 5% obtained when the calorimeters 
alone are used for jet clustering. 

Neutrinos escape the detector undetected and give rise to the missing transverse momentum 
vector, defined as the projection on the plane perpendicular to the beams of the negative vector 
sum of the momenta of all reconstructed particles in an event. Its magnitude is referred to as 

^miss 

The jets contain neutral particles from pileup events. The contribution from these additional 
particles is subtracted based on the average expectation of the energy deposited from pileup in 
the jet area, using the methods described in Ref. Il43l . 

Eor the identification of jets resulting from fragmentation of b quarks ("b jets"), an algorithm is 
used that combines information from reconstructed tracks and from secondary vertices, both 
caracterized by a displacement with respect to the PV. This information is combined into a sin¬ 
gle discriminating variable and jets are tagged as b jets based on its value. The algorithm is 
referred to as "combined secondary vertex tagger" and is described in Ref. Il44l . In most of 
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the analyses described in the following, a minimum value of this variable (medium operating 
point) is chosen such that the b tagging efficiency is 70% and th e lig ht-flavor jet misidentifica- 
tion rate is 1% in tf events. The analyses presented in Sections 5.2 and 5.5 also use a smaller 
minimum value of the discriminating variable (loose operating point), yielding a higher effi¬ 
ciency of approximately 80%, with a light-flavor misidentification rate of 10%. 


4.1 Jet substructure methods 

Because of the possible large mass of the T quarks, the top quarks, Higgs and Wbosons from 
T quark decays might have significant Lorentz boosts. Daughter particles produced in these 
decays would therefore not be well separated. In many cases, all decay products are clustered 
into a single large jet by the event reconstruction algorithms. These merged jets exhibit an 
intrinsic substructure that can be analyzed with dedicated jet substructure algorithms. In order 
to cluster the decay products from top quarks and Higgs boson into wide jets, the CA algorithm 
is used with size parameters R=1.5 (CA15 jets) or R=0.8 (CA8 jets). A number of jet substructure 
algorithms are then used in different analyses to identify jets from top quark or Higgs boson 
decays. This process is known as t or H tagging, and in some cases relies on b tagging of 
individual subjets. 

The inclusive T quark search in final states with leptons discussed in Section [5^ uses the CM- 
STopTagger ||45]| , which is based on the algorithm developed in Ref. Eb]. The tagger identi¬ 
fies a top quark decay if a CA8 jet with pj > 400 GeV is found with a mass between 140 and 
250 GeV and at least three subjets with a minimum mass of subjet pairs larger than 50 GeV. The 
sensitivity of the GMSTopTagger is suitable for a regime with jet pj > 400 GeV where the 
decay products are collimated to be within the acceptance of a jet with the size parameter of 
0 . 8 . 

The search for T —)■ tH in the hadronic final state (Section |5.3[ ) adopts the HEPTopTagger al¬ 
gorithm Il47ll48l . which employs GA15 jets to increase the acceptance to top quarks with a mod¬ 
erate Lorentz boost (pj > 200 GeV). This facilitates a smooth transition between the boosted 
and resolved regimes. A GA15 t jet candidate is required to exhibit a substructure compatible 
with a three-body decay. If this requirement is satisfied, the HEPTopTagger clustering algo¬ 
rithm identifies the three subjets, and then requires that the mass of a subjet pair be consistent 
with the Wboson mass and the mass of the three subjets be consistent with the top mass. The 
t tagging performance is further enhanced by the application of b tagging to subjets of GA15 
jets [49|. Subjet b tagging is also used to identify decays of boosted Higgs bosons into a bottom 
quark-antiquark pair. The subjets of GA15 jets are reconstructed using the filtering algorithm 
described in Ref. Il50l . Two filtered subjets of GA15 jets are required to have a di-subjet invariant 
mass larger than 60 GeV. Both subjets are tagged using the subjet b tagging algorithm, which is 
based on the same algorithm used for regular anti-kj jets, discussed above, with the difference 
that only tracks and secondary vertices associated with the individual subjets are used to build 
the b tag discriminator. 

Por the identification of boosted Wbosons, two subjets are required to be reconstructed by a 
pruning algorithm |[50lj52l . The mass of the pruned jet has to be compatible with the mass of 
the Wboson, within a mass window that differs slightly depending on the analysis considered. 
The inclusive analysis in Section [5d] requires a Wjet to have pj > 200 GeV and a mass between 
60 and 130 GeV. The search for T —> bW with single leptons (Section |5.2[ ) applies the same 
Pj selection, but the mass window is tightened to 60 to 100 GeV. The search for T —> bW in 
hadronic final states (Section |5.4[ | requires pj > 150 GeV in combination with a jet mass mj 
requirement of 60 < nij < 100 GeV. Additionally, this analysis complements pruning with a 
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5 Analysis channels 


selection on the mass drop 1501 , which is defined as the ratio of the largest subjet mass to that 
of the original jet. Requiring the mass drop to be <0.4 rejects events containing massive jets 
from QCD multijet processes. 

The different performance of the t tagging and Wtagging algorithms in data and simulation is 
taken into account with scale factors that are applied to the simulated events Il48ll53l . 


5 Analysis channels 

In this Section, five distinct searches for T quarks are presented, each optimized for a different 
topology. The analyses described in Sections 5.1 and 5.2 are based on leptonic final states. 
While the former is an inclusive search covering all possible decay modes, the latter is a search 
specifically optimized to find T —> bW decays. The searches presented in Section 5.3 and 
Section 5.4 are optimized for boosted event topologies in hadronic final states and make use 
of jet substructure techniques. Finally, the search treated in Section 5.5 is sensitive to T ^ tH 
decays, where the Higgs boson decays to a pair of photons. 


5.1 Inclusive search with singie and muitipie ieptons 

The inclusive search described in this Section is sensitive to all decay modes of the T quark, i.e., 
T —> tH, T —> tZ, and T —> bW. It is divided into two channels: one charmel in which exactly 
one lepton is selected and the other channel with at least two leptons. Further details are given 
in Ref. ITFI . 


5.1.1 Single-lepton channel 

Single-lepton events must contain exactly one isolated muon or electron with pj > 32 GeV. 
In addition to the lepton, events must also have at least three AK5 jets with pj > 120, 90, 
and 50 GeV. A fourth AK5 jet with pj > 35 GeV is required if no Wjet is identified in the event. 
To fulfill the lepton isolation requirement, jets must be separated by AR > 0.4 from muons and 
by AR > 0.3 from electrons. The requirement on the jet multiplicity and pj significantly sup¬ 
presses background processes. The contribution from QGD multijet events is further reduced 
by selecting events with > 20 GeV. The major selection requirements are summarized in 
Table |2] 

Table 2: Main selection requirements for the single-lepton analysis. 


Variable 

Selection 

Pt lepton 

>32 GeV 

Number of jets 

>3 

Pt jets 

>120, 90, and 50 GeV 

Wtag 

>1 or >1 jets with pj > 35 GeV 

^miss 

>20 GeV 


Some background events from W+jets production remain after the event selection. This process 
is not well modeled by simulations and the normalization is determined from a control sample 
in data. This sample is defined by single-lepton events fulfilling the signal selection criteria, but 
failing the requirement that a fourth jet with pj > 35 GeV or alternatively a Wjet is identified 
in the event. 

A boosted decision tree (BDT) f54l is used to discriminate between signal and background 
events. Different BDTs are implemented for events with and without identified W jets and for 
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each hypothetical value of the mass of the T quark. The use of dedicated BDTs for different T 
quark decay modes does not improve the performance, so the BDTs are trained irrespective of 
the branching fraction of the T quark. 

The variables used for the calculation of the BDT discriminant are jet multiplicity, b-tagged jet 
multiplicity, lepton pj, pj of the third jet, pj of the fourth jet, and Hj, where Hj is de¬ 

fined as the scalar pj sum of all jets with pj > 30 GeV. For events with at least one Wjet, the 
multiplicity and pj of W-tagged jets and the numbers of t-tagged jets are also included in the 
BDT training. These variables are chosen based on their discrimination power as calculated by 
the BDT algorithm, and on the absence of significant correlations between the different vari¬ 
ables. The final BDT distributions are shown in Ref. flSH . The total numbers of events predicted 
for background processes and observed in collision data are shown in Table The predicted 
contributions for each background process are available in Ref. fl^ . The signal selection effi¬ 
ciencies are between 7.5% and 9.4% which corresponds to an expected number of 850 events 
for a T quark mass of 500 GeV and 6 events for a T quark mass of 1000 GeV assuming branching 
fractions to tH, tZ, and bW of 25%, 25%, and 50%, respectively A detailed table with selection 
efficiencies and expected number of events is available in Ref. flSl . 

Table 3: Numbers of events predicted for background processes and observed in collision data 
for the single-lepton analysis. The uncertainties include those in the luminosity, the cross sec¬ 
tions and the correction factors on lepton and trigger efficiencies. From Ref. KTSlI . 



Muon 

Electron 

Total background 

61900 ± 13900 

61500 ± 13700 

Data 

58478 

57743 


5.1.2 Multilepton channel 

This channel uses four mutually exclusive subsamples with at least two leptons: two opposite- 
sign dilepton samples (referred to as OSl and OS2 samples) which differ by the required num¬ 
bers of jets in the event, a same-sign dilepton sample (the SS sample) and a multilepton sample. 
The division into opposite- and same-sign dilepton events is based on the charge of the leptons. 

Multilepton events must contain at least three leptons with pj > 20 GeV. To reject backgrounds 
from heavy-flavor resonances and low-mass Drell-Yan (DY) production, multilepton events 
must contain a dilepton pair of the same flavor and of opposite charge with an invariant mass 
above 20 GeV. Events in which < 30 GeV are discarded. Jets must be separated by AR > 
0.3 from the selected leptons and at least one of the jets has to fulfill the b tagging criteria. 

The OSl dilepton sample targets events in which both T quarks decay to bW fl^ . This dilepton 
sample contains events with either two or three jets, Hj > 300 GeV, and Sj > 900 GeV, where 
St is the sum of Hj, £^^®®, and the transverse momenta of all leptons. Events are discarded 
where there is a dilepton pair with same-flavor leptons and a mass consistent with that 
of a Zboson (76 < Mu < 106GeV). To reduce the tt background, all the possible pair-wise 
combinations of a lepton and a b jet are considered and their invariant masses are all required 
to be larger than 170 GeV. 

The DY background is not modeled reliably in the selected kinematic region and is controlled 
using a data sample consisting of events with no b-tagged jets, £3^®® <10 GeV, Sj < 700 GeV, 
and Ht > 300 GeV. 

The OS2 dilepton sample consists of events with at least five jets, two of which must be identi¬ 
fied as b jets. Events are also required to have Hj > 500 GeV, and Sj > 1000 GeV. This sample 
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is mostly sensitive to signal events where both T quarks decay to tZ. The dominant background 
is tt production. 

The SS sample selection criteria target events in which at least one T quark decays to tZ or 
tH. Besides the lepton selection criteria, at least three jets are required, Hj > 500 GeV, and 
St > 700 GeV. 

Different processes contribute to the background in the SS sample. A minor contribution is 
given by SM processes leading to prompt SS dilepton signatures, which have very small cross 
sections. These processes can be simulated reliably. The prompt OS dilepton production can 
also contribute if one lepton is misreconstructed with the wrong sign of the charge. The misre- 
construction probability of the charge sign is negligible for muons in the kinematic range con¬ 
sidered, while for electrons it is determined from control data samples. We determine the prob¬ 
ability to misreconstruct the charge sign of an electron from events with a dileptonic Zdecay, 
selected with the same criteria as in the signal selection except for the charge requirement. In¬ 
strumental backgrounds in which misidentified jets create lepton candidates are determined 
from control data samples in which non-prompt and fake leptons are enriched. 

The multilepton sample, like the SS sample, is mostly sensitive to signal events in which at 
least one T quark decays to tZ or tH. The backgrounds are suppressed by selecting events with 
at least three jets, Hj > 500 GeV, and Sj > 700 GeV. Prompt backgrounds in this channel are 
due to SM processes with three or more leptons in the final state, such as diboson and triboson 
production. These are correctly modeled by simulation. Nonprompt backgrounds are caused 
by the misidentification of one or more leptons, by tt production, and by other processes. As for 
the dilepton samples, data control samples are used to evaluate these sources of background. 

The main selection requirements for the four samples are summarized in Table 

Table 4: Main selection requirements for the four multilepton channels: the opposite-sign dilep¬ 
ton samples with two or three jets (OSl) and with at least five jets (OS2), the same-sign dilep¬ 
ton sample (SS), and the multilepton sample. The smallest mass obtained from all the possible 
combinations of leptons and b jets is indicated by Mb^. 



OSl 

OS2 

SS 

Multileptons 

Ht (GeV) 

>300 

>500 

>500 

>500 

St (GeV) 

>900 

>1000 

>700 

>700 

Number of jets 

2 or 3 

>5 

>3 

>3 

b tags 

>1 

>2 

>1 

>1 

(GeV) 

>30 

>30 

>30 

>30 

Mm (GeV) 

>170 

— 

— 

— 

(GeV) 

>20 

>20 

>20 

>20 

Z boson veto 

yes 

no 

no 

no 


The numbers of events in the multilepton samples are given in Table both for data and for 
estimated background contributions. The predicted contributions for each background pro¬ 
cess are available in Ref. IfTSl . The selection efficiencies for signal events are between 0.15% 
and 0.44% which corresponds to an expected number of 16.7 events for a T quark mass of 
500 GeV and 0.28 events for a T quark mass of 1000 GeV, assuming branching fractions to tH, 
tZ, and bW of 25%, 25%, and 50%, respectively. A detailed table with selection efficiencies and 
expected number of events is available in Ref. IITSl . The numbers of background and signal 
events are of similar order of magnitude. The sensitivity to the signal is enhanced by further 
splitting the samples according to the lepton flavor. The dilepton samples are separated into 
three subsamples, fifi, f/e, and ee. The multilepton sample is divided into a subsample. 
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an eee subsample, and a third subsample with events with mixed lepton flavors. Data and SM 
background expectations are found to be in agreement. 

Table 5: Numbers of events selected in data and expected for the backgrounds. Shown are 
the opposite-sign dilepton samples with two or three jets (OSl) and with at least 5 jets (OS2), 
the same-sign dilepton sample (SS), and the multilepton sample. The background sources not 
contributing to the channel are indicated by a dash ("-")■ The uncertainties include statistical, 
normalization, and luminosity uncertainties. From Ref. Ill5l . 



OSl 

OS2 

SS 

Multileptons 

Total background 

17.4 ± 3.7 

84 ± 12 

16.5 ± 4.8 

3.7 ± 1.3 

Data 

20 

86 

18 

2 


5.2 Search for T bW with single leptons 

The analysis described in this Section is optimized for the event topology in which both T 
quarks decay into a bottom quark and a Wboson. 

Events are required to have one isolated muon or electron, where muon candidates must have 
Pj > 45 GeV and electron candidates must have pj > 30 GeV. At least four jets are required, 
either at least four AK5 jets or at least three AK5 jets plus at least one GAS jet. The AK5 jets are 
required to have pj > 30 GeV and GAS jets are required to have pj > 200 GeV. Both types of 
jets must have \fl\ < 2.4. 

The GAS jets are used to identify merged hadronic decays of Wbosons with high Lorentz boost. 
The AK5 jets are replaced by the two pruned subjets of W-tagged GAS jets if the angular dis¬ 
tance between AK5 and GAS jets fulfills the matching criterion AR(JetcA8/ JetAKs) < 0.04. Un¬ 
matched AK5 jets and the subjets of matched W-tagged GAS jets are used as input for a kine¬ 
matic fit, which is described below. The four jets or subjets are required to satisfy pj > 120,90, 
50, and 30 GeV. At least one of the AK5 jets has to satisfy the b tagging criteria. 

A kinematic fit is made to each event for the hypothesis TT —> bW+bW“ —)■ f'vbqqh, subject 
to the constraints, m{lv) = wj(qq') = Mw/ and m{ivb) = ?«(qq'b) = Mfn, the fitted mass of 
the selected T candidate. The in the event is attributed to the undetected neutrino from 
leptonic Wdecays. If a selected event has more than four jets, the fifth jet with highest px is also 
considered and all the possible combinations of four jets are tested in the kinematic fit. 

Only events containing fit combinations with probability p{x^) > 1% are retained. The 
efficiency of the p{x^) criterion is 62% for signal events with a T quark mass of 800 GeV while 
76% of background events are rejected. The p{x^) criterion removes badly reconstructed events 
with poor mass resolution and improves the signal-to-background ratio in the reconstructed 
mass spectrum. 

To reduce the large combinatorial background, the b tagging and the Wtagging information is 
used. If a Wtag is present, only those combinations where the subjets of the Wjet match the 
Wdecay products are considered. The best combination is selected from groups of fit combina¬ 
tions with decreasing b tag multiplicity, ranked by the b tagging operating point (OP), as listed 
below: 

• 2b tags at medium OP; 

• lb tag at medium OP and 1 b tag at loose OP; 

• lb tag at medium OP; 

• 2b tags at loose OP. 
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5 Analysis channels 


Decay products of T quarks have on average higher pj than those from the SM backgrounds. 
To suppress the backgrounds and enhance the signal significance, we select events with large 
values of the St variable, which is defined here as a sum of pj of the lepton, and pj of the 
four jets that minimize the io the kinematic fit. Figure [^demonstrates that SM backgrounds 
and a T quark signal populate different regions in the two-dimensional Sj and Mgt distribution. 
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Figure 1: Correlation between the Sj and the M^t observables in the search for T —> bW with 
single leptons, for background processes (left) and for a simulated signal, with a T quark mass 
of 800 GeV (right). The color gradient indicates the entries per bin in arbitrary units (a.u.). 


We test the modelling of the shape of the reconstructed mass, and verify how well the SM 
background expectations agree with data, as a function of Sj. Figure j^shows the reconstructed 
mass distributions separately for p+jets and e+jets events with the Sj > 1000 GeV requirement. 
Gorrectly reconstructed tt events peak near the top quark mass value, while events with mis- 
assigned jets constitute a combinatorial background, and populate a region of higher masses, 
where the potential signal is expected to appear. Table [^ (left columns) presents the event 
yields of SM backgrounds and data for this selection. The dominant background process is tt 
production. Smaller but still significant backgrounds come from W+jets and single top quark 
production. In the e+jets channel there is also a contribution from QGD multijet production. 
Other backgrounds have been found to be negligible. Data and SM background expectations 
agree in both shape and total normalization. 

Table 6: Numbers of observed and expected background events after the event selection. The 
uncertainties in the predicted numbers of events include both the statistical and systematic 
uncertainties. 



Selection (Sj > 1000 GeV) 

Selection (St > 1240 GeV) 


p+jets 

e+jets 

p+jets 

e+jets 

tt 

325 ± 37 

279 ± 35 

51 ±6 

52 ±6 

W + >3jets 

49 ±8 

60 ±9 

18 ±3 

19 ±4 

Single top 

20 ±5 

36 ±10 

6.9 ±2.3 

10 ±4 

Z/ 7 * -|- >3jets 

3.9 ±0.8 

3.3 ±0.6 

1.4 ±0.4 

1.1 ±0.3 

WW, WZ, ZZ 

3.1 ±1.0 

<1 

<1 

<1 

Multijet 

<1 

18 ±4 

<1 

6.1 ±1.7 

Total background 

401 ± 38 

396 ± 38 

77 ±7 

88 ±9 

Data 

417 

398 

81 

83 


We apply a requirement of Sj > 1240 GeV in the final event selection. This condition is opti¬ 
mized to enhance the sensitivity to the signal, based on SM backgrounds and T signal expecta- 











5.3 All-hadronic search for T tH 


11 



Figure 2: Distributions of the reconstructed mass M^t for /r-i-jets (left) and e-i-jets (right) events. 
The data are shown as points and the simulated backgrounds as shaded histograms. The 
hatched region and the shaded area in the lower panel represent the statistical uncertainty 
in the background. The expected signal (dotted line) for a T quark with a mass of 800 GeV is 
multiplied by a factor of 5 for better visibility. The lower panel represents the ratio between 
data and the sum of the backgrounds (BG). The overflow of the distributions is added to the 
last bin. 

Table 7: Main selection requirements for the T —> bW search with single leptons. 


Variable 

Selection 

Pj muon 

>45 GeV 

Pt electron 

>30 GeV 

Number of jets 

>4 

Pt jets 

>120, 90, 50, and 30 GeV 

Wtags 

0 or 1 

b tags 

1 or 2 

St 

>1240 GeV 

^miss 

>30 GeV 


tions. The major selection requirements are summarized in Table 

Table|^(right columns) presents the event yields for expected SM backgrounds and data. Signal 
efficiencies are of the order of 0.5-4% for T quark masses from 500 to 1000 GeV. They are 
summarized in Table |8l 

The Mfit distribution for the final event selection is shown in Fig. The f/-i-jets and e-i-jets final 
states give very similar results. The observed data are compatible with background expecta¬ 
tions from SM processes. The /r-i-jets and e-i-jets channels are combined to improve the statistics 
for the simulated SM backgrounds. 

5.3 All-hadronic search for T —tH 

This channel is optimized for the event topology in which at least one T quark decays to T —> 
tH, where the top quark decays into bW and the Wboson decays hadronically, and the Higgs 
boson decays into two b quarks. Because of the expected high mass of the T quarks, the top 
quarks and Higgs bosons can have significant Lorentz boost; therefore the event selection is 
based on jet substructure requirements, as described in Section [hi] 

At least one t-tagged and one H-tagged GA15 jet are required, where the t-tagged jets must 
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5 Analysis channels 


Table 8: Selection efficiencies and numbers of expected signal events for the selection Sj > 
1240 GeV, for the two channels of the T —> bW search with single leptons. Different T quark 
mass hypotheses are considered and a 100% branching fraction to bW is assumed. 

T quark mass Muon channel Electron channel 
(GeV) Efficiency Events Efficiency Events 
500 
600 
700 
800 
900 
1000 
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Eigure 3: Distributions of the reconstructed T quark mass Mfit for bWhW candidate events in 
the search for T —)■ bW with single leptons, combining the /r+jets and e+jets samples after the 
selection Sj > 1240 GeV. Data are shown as points and the simulated backgrounds as shaded 
histograms. The hatched region and the shaded area in the lower panel represent both the 
statistical and the systematic uncertainties in the total background. The expected signal for a 
T quark of mass 800 GeV is multiplied by a factor of 2. The lower panel represents the ratio 
between data and the sum of the backgrounds (BG). The horizontal error bars represent the bin 
width. The overflow of the distribution is added to the last bin. 

have pj > 200 GeV and the H-tagged jets must have pj > 150 GeV. Two variables are used 
to further distinguish the signal from the background events after the event selection. These 
variables are defined here as the scalar pj sum of subjets of GA15 jets, and the invariant 
mass of two b-tagged subjets in the H-tagged jets. These two variables are used for setting 
upper limits on the T quark production cross section. The major selection requirements are 
summarized in Table HI 

Backgrounds due to QGD multijet production are determined from data using signal-depleted 
sideband regions. These sidebands are defined by inverting the jet substructure criteria. Back¬ 
grounds due to tf events are determined from simulation; other backgrounds are found to be 
negligible. 

To maximize the sensitivity of the analysis, the events are divided into two categories: a cat¬ 
egory with a single H tag and a category with at least two H tags. The background estimates 
are well matched to the observed data, as discussed in Ref. HU. Eor the final event selection, 
the and variables are combined into a single discriminator using a likelihood ratio 
method. The numbers of expected background events and events observed in data after the 


0.50% 

59 

0.46% 

53 

1.24% 

43 

1.30% 

44 

2.38% 

28 

2.38% 

27 

3.04% 

13 

3.17% 

13 

3.48% 

5.6 

3.63% 

5.8 

3.52% 

2.3 

3.86% 

2.5 
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Table 9: Main selection requirements for the all-hadronic search for T —> tH. 


Variable 

Selection 

LTSUb 

>720 GeV 

Number of GA15 jets 

>2 

Pt GA15 jets 

>150 GeV 

Pt t-tagged jets 

>200 GeV 

Number of t tags 

>1 

Number of H tags 

>1 


full selection are shown in Table [fO] The observed data are compatible with background expec¬ 
tations from SM processes. The signal selection efficiencies are between 2.5% and 7.2% which 
corresponds to an expected number of 283 signal events for a T quark mass of 500 GeV and 4.9 
events for a T quark mass of 1000 GeV, assuming ,B(T —> tH) = 100%. A detailed table with 
selection efficiencies and expected numbers of signal events is available in Ref. IIT9I . 

Table 10: Predicted numbers of total background events and observed events for the two event 
categories with one and with multiple H tags, for the all-hadronic search for T —> tH. The 
quoted uncertainties are statistical only. From Ref. Ill9]| . 



Single H tag category 

Multiple H tags category 

Total background 

1403 ± 14 

182 ± 5 

Data 

1355 

205 


5.4 All-hadronic search for T bW 

This channel is optimized for the event topology in which both T quarks decay to T —> bW, 
where the W bosons decay hadronically Events are selected by requiring two W-tagged GAS 
jets with pj > 150 GeV. At least two additional AK5 jets with pj > 50 GeV are required, one of 
which must be b-tagged. Events are divided into categories defined by the numbers of b-tagged 
jets: one or at least two. 

After the event selection, two T candidates Ti and T 2 are reconstructed using combinations of 
the Wjets and the AK5 jets. The order of Ti and T 2 is arbitrary. The reconstruction is performed 
by identifying the combinations of Wjets and AK5 jets having the smallest invariant mass dif¬ 
ference. Eigure|^ shows the two-dimensional distribution of the masses of each reconstructed 
T candidate in a signal sample with a simulated T quark mass of 800 GeV. The reconstructed 
mass peak is clearly visible at the expected value. The misreconstruction rate, where the wrong 
combination of jets is chosen, is small and does not affect the signal acceptance. Additional 
event requirements are then applied to increase sensitivity to the signal process. The T can¬ 
didate masses must be greater than 200 GeV, and the fractional difference Uf in the masses of 
the two T candidates m(Ti) and m{T 2 ), where Uf = |m(Ti) — m{T 2 )\/ (m(Ti) -|- m(T 2 )), must 
be less than 10%. The two T candidates must fall in opposite hemispheres of the detector, 
A(^(Ti,T 2 ) > 5n/6, and finally must be above 1000 GeV, where is defined as the 
scalar pj sum of the four jets used to reconstruct the T candidates. The major selection require¬ 
ments are summarized in Table [TT] 

The dominant backgrounds are due to QGD multijet production and tt production. Other back¬ 
ground contributions are negligible. 

To obtain the shape of the QGD multijet background, a control region is defined by requir¬ 
ing > 1000 GeV, but inverting the requirement on the fractional mass difference, > 0.1. 
This control region is enriched in multijet events and has a negligible signal contamination. The 
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5 Analysis channels 



Mass of the candidate (GeV) 

Figure 4: Two-dimensional distribution of the masses of each reconstructed T candidate in the 
selected events for the all-hadronic search for T —> bW, for a simulated signal sample with a T 
quark mass of 800 GeV. The order of Ti and T 2 is arbitrary. 

Table 11: Main selection requirements for the all-hadronic search for T —)■ bW. 


Variable 

Selection 

Number of AK5 jets 

>2 

pj AK5 jets 

>50 GeV 

Number of W-tagged jets 

>2 

Pj W-tagged jets 

>150 GeV 

Reconstructed T candidate mass 

>200 GeV 

a/ 

<10% 

A(/>(Ti,T2) 

>571/6 

1x4 jet 

>1000GeV 


shape of the distribution in the control region, after subtracting the expected tt contribu¬ 
tion, is used to model the QCD multijet events entering the signal region. The distribution 
in the signal region agrees with the distribution in the sideband region for simulated QCD mul¬ 
tijet events. The normalization of the QCD multijet background is not fixed, and is determined 
in the limit setting procedure. This procedure is done independently for events containing one 
and at least two b-tagged jets. 


Figure]^ shows the post-fit H^extAjet distributions obtained with the above method. Data are 
found to be in agreement with the expected background contributions. The numbers of ex¬ 
pected background events and events observed in data after full selection are shown in Ta¬ 
ble The numbers of expected signal events and selection efficiencies assuming i8(T 


bW) = 100% are summarized in Table 13 


5.5 Search for T —tH with H —77 

The analysis described in this section is optimized for events with one T quark decaying to tH, 
where the Higgs boson decays into a pair of photons. The main advantage of this channel is 
the possibility to precisely measure the invariant mass of the diphoton system (niyy) so that a 
peak in the nijj distribution would be present for signal events. The disadvantage is the small 
Higgs branching fraction of the order of 2 x 10^^ liSSl . The analysis concept is the same as for 
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Figure 5: The distributions for single b tag events (left) and for events with at least two 
b tags (right) for the all-hadronic search for T —)■ bW, including the QCD multijet background 
estimate obtained from data and the T quark signal with a mass of 800 GeV. The hatched 
region and the shaded area in the lower panel represent both the statistical and the systematic 
uncertainties in the total background. The lower panel represents the ratio between data and 
the sum of the backgrounds (BG). The horizontal error bars represent the bin width. 


Table 12: Summary of expected and observed background yields for the two channels of the 
T —> bW search in the all-hadronic final state. 



1 b tag channel 

>2 b tags channel 

tt 

20.3 ± 1.3 

3.45 ± 0.55 

QGD multijet 

979 ± 29 

80.2 ± 6.4 

Total background 

999 ± 31 

84 ±7 

Data 

998 

84 


searches of the SM Higgs boson in the H ^ 77 decay channel Il56l . 

Events with two isolated photons are selected. Additional leptons and jets coming from the 
decay of top quarks or a second Higgs boson are required. In order to maximize the sensitivity 
of the analysis, two search channels are defined, targeting different decay modes of the top 
quark: 

• the leptonic channel searches for events with a pair of photons and at least one iso¬ 
lated high-px muon or electron; 

• the hadronic channel searches for events with a pair of photons and no isolated 
muons or electrons. 

The resonant contributions from the tfH background are determined from simulation. The 
nonresonant contribution is composed of events with two prompt photons arising from QGD 
multijet production as well as for emission in top quark production ( 77 +jets, tf -|- 77 , t -|- 77 ). 
The tt events are more likely to have a jet misreconstructed as a photon, because of the large 
numbers of jets in the final state. The simulation of such sources of instrumental background 
is not completely reliable. The background model is therefore derived from data. 

The control sample used to estimate the nonresonant background consists of events where at 
least one photon passes loose identification requirements but does not pass the final event selec¬ 
tion. This sample is enriched with events containing one misidentified photon. A reweighting 
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6 Combination strategy 


Table 13: Selection efficiencies and numbers of expected signal events, for the two channels 
of the T —> bW search in the hadronic final state. Different T quark mass hypotheses are 
considered and a 100% branching fraction to bW is assumed. 


T quark mass 
(GeV) 

1 b tag channel 
Efficiency Events 

>2 b tags channel 
Efficiency Events 

500 

1 . 01 % 

103.4 

0 . 86 % 

84.7 

600 

2.24% 

66.0 

1.81% 

52.5 

700 

3.15% 

31.24 

2.35% 

22.80 

800 

4.07% 

14.64 

2.51% 

8.79 

900 

4.68% 

6.57 

2.44% 

3.33 

1000 

4.95% 

2.81 

2.35% 

1.29 


is applied, in order to match the px and p spectra of the photons in this control sample to those 
obtained after the signal selection. This is done independently for each photon. 

The event selection is based upon six quantities that have the largest discriminating powers 
between signal and backgrounds and that have small correlations. They include the transverse 
momenta of the larger px photon ( 71 ) and smaller px photon ( 72 ). The selection criteria are 
optimized to produce the most stringent limits on the signal cross section and are listed Table [l^ 
for both leptonic and hadronic channels. 

Table 14: Final selection criteria for hadronic and leptonic channels of the search for T —> tH 
with H —> 77 . 


Variable 

Leptonic channel 

Hadronic channel 

Et(7i) 

>\myy 

>imyy 

Prill) 

25 GeV 

35 GeV 

Number of jets 

>2 

>2 

Sx 

>770 GeV 

>1000GeV 

Leptons 

>1 

0 

b tags 

— 

>1 


The nonresonant background contributions are obtained from unbinned maximum likelihood 
fits to the diphoton mass distribution over the range 100 < < 180 GeV, under the hypothe¬ 

sis of no signal. An exponential function is chosen for these fits. Studies of pseudo-experiments 
showed that the use of an exponential function does not introduce a bias in the estimation of the 
numbers of background events in both categories. In Fig. the observed diphoton mass dis¬ 
tribution in each event category is shown, together with the expected signal and the expected 
resonant background contribution. The error bands show the uncertainty in the background 
shapes associated with the statistical uncertainties of the fits. The numbers of expected back¬ 
ground events and events observed in data after final selection are shown in Table The 
numbers of expected signal events and selection efficiencies assuming S(T —)■ tH) = 100% are 
summarized in Tabled^ 

The data in the signal window are compatible with background expectations from SM pro¬ 
cesses. 


6 Combination strategy 

The event samples selected by the five analyses are almost entirely distinct and therefore, sig¬ 
nal limits extracted from those analyses are statistically independent. They can be combined to 
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Figure 6 : Diphoton invariant mass distribution for the leptonic (left) and hadronic (right) chan¬ 
nels of the search for T —> tH with H—> 77 . The signal is normalized to the predicted theoretical 
cross section corresponding to mj = 700 GeV. The backgrounds predicted by the fit are shown 
as a solid line while the corresponding uncertainties are shown as bands around the line, where 
the inner band indicates the In and the outer band indicates the 2n uncertainties. Bins with zero 
entries are not shown. 


Table 15: Expected yields for ttH and nonresonant background (from the fit to data) and the 
numbers of observed events in data after full event selection for the two channels of the T —)■ tH 
search in the final state with photons. All the yields are computed in a window of 1 full width 
at half maximum i.e., 125 ±1.5 GeV. 



Leptonic channel 

Hadronic channel 

ttH 

U.UDV_QQgg 

U.U^z_o.o06 

Nonresonant background 

n 11 +0.07 
^•^'*^-0.03 

0 

^•"^-0.13 

Total background 

Q 1 c+0.07 

^•"^-0.13 

Data 

0 

2 


yield a result that is more stringent than any of the inputs. Because the backgrounds are largely 
common to all analyses, the background estimates are largely correlated but well determined 
by the multiple independent samples. In particular, most analyses have top quark pair produc¬ 
tion as a background process. This background normalization is correlated among the analyses 
in the combination, providing for the combination a better background estimation than in the 
individual analyses. Similar arguments hold for the correlated systematic uncertainties, which 
are discussed in more detail in Section l6d1 


The inclusive analysis with single and multiple leptons described in Section 5.1 is able to set 
limits for all T quark decay modes. Dedicated optimizations to enhance the sensitivity for 
T ^ bW decays are described in Section [5^ These optimizations use single-lepton events. To 
avoid double counting of events we replace the single-lepton part of the inclusive approach 
(Section |5.1 1 with the single-lepton analysis described in Section [5^ This is done for scenarios 
with jB(T —> bW) values of at least 80%. For lower S(T —>• bW) values this approach is inferior 
and we use the inclusive results from Section [5d] only. At every point the approach used is that 
which gives the best expected limit. The other three analyses described in Sections [53| to |5.5| do 
not have any overlap so they are always combined with the cases above. 
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6 Combination strategy 


Table 16: Selection efficiencies and numbers of expected signal events, for the two channels 
of the T —> tH search in the final state with photons. Different T quark mass hypotheses are 
considered and a 100% branching fraction to tH is assumed. 


T quark mass 
(GeV) 

Leptonic channel 
Efficiency Events 

Hadronic channel 
Efficiency Events 

500 

6.7% 

6.0 

9.3% 

8.3 

600 

9.6% 

8.7 

18.1% 

16.4 

700 

11 . 0 % 

9.8 

26.0% 

23.8 

800 

12 . 0 % 

10.9 

30.0% 

27.3 

900 

11.4% 

10.4 

32.0% 

29.3 


For the statistical combination a Bayesian method 157 1 has been adopted in which the sys¬ 
tematic uncertainties are taken into account as nuisance parameters with their corresponding 
priors as discussed in Section 6.1 Upper limits on the T quark production cross section are 


obtained with the Theta framework f58 l. Systematic uncertainties are taken into account as 
global normali z ation uncertainties and as shape uncertainties where applicable. More details 
about the treatment of systematic uncertainties are given in the next section. 


6.1 Systematic uncertainties 

Some of the individual analyses are sensitive to the same systematic uncertainties, for example 
the uncertainty in the integrated luminosity, the jet energy scale and the b tagging efficiency. 
Such uncertainties are treated as fully correlated, as is done technically by correlating the corre¬ 
sponding nuisance parameters in the limit setting procedure. This treatment allows improved 
constraints to be obtained on these parameters than is possible in the standard analyses. 

The systematic uncertainties fall into two types: those which affect the normalization of the 
signal and background samples, and those which also affect the shapes of distributions. The 
uncertainty in the tf cross section is 13%. It is obtained from the tf cross section measure¬ 
ment for large invariant mass values of the tf system. The uncertainty in the integrated 
luminosity is 2.6% Il60| . 

Shape uncertainties include the jet energy scale, the jet energy resolution and the b tagging 
efficiency uncertainties. We also consider the uncertainties in the efficiencies of the t tagging, 
Wtagging, and H tagging algorithms Il4^l4^l53l . The uncertainty due to the energy deposits 
not associated with jets (unclustered energy) has an impact on the missing pj. This effect is 
taken into account in the single-lepton channel. The size of this uncertainty typically varies 
from a few percent up to 10 %. 

The systematic uncertainty in the pileup jet identification is taken into account in the analysis 
with H —> 77 . It is derived through the use of the data/simulation scale factors (SF), which are 
binned in jet rj and pj Il56l . 

For the photon identification efficiency, the uncertainty in the SF is taken into account. The 
SF corrects the efficiency in simulation to the efficiency as measured in data using a "tag-and- 
probe" technique 1611 applied to Z —> e+e^ events. The uncertainty applied to this SF amounts 
to 3% in the barrel region of the calorimeter and 4% in the endcaps. 

Lepton trigger efficiencies, lepton identification efficiencies, and corresponding correction fac¬ 
tors for simulated events are obtained from data using decays of Zbosons to dileptons. These 
uncertainties are < 3%. 

For simulated tf and ttH events, uncertainties due to renormalization and factorization scales 
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(f/R and ^f) are taken into account by varying both scales simultaneously up and down by a 
factor of two. Uncertainties arising from the choice of PDFs are taken into account. Simulated 
background events are weighted according to the uncertainties parameterized by the CTEQ6 
eigenvectors ISTI . The shifts produced by the individual eigenvectors are added in quadrature 
in each bin of the relevant distributions. 


A systematic uncertainty of 50% is assigned to the diboson backgrounds, single top quark 
production and the Wand Zboson background. This accounts for the effects of the /^r and fip 
variations in simulation and the uncertainties in the determination of the W+jets SF from data. 


Modified "template" distributions of those quantities that are affected by the respective un¬ 
certainties are obtained by varying the respective quantity by its uncertainty, namely by ±1 
standard deviation. In the limit setting procedure a likelihood fit is performed in which the 
nominal distribution and the modified templates are interpolated. The corresponding uncer¬ 
tainty is represented as a nuisance parameter, which receives its prior constraints from the tem¬ 
plate distributions. In the fit, the templates are allowed to be extrapolated beyond ±1 standard 
deviation, but this happens rarely The resulting fit values are always within ±1.5 standard 
deviations of their prior values. 


The list of nuisance parameters of all analysis channels is shown in Table 17 This table also 
indicates which parameters are correlated and which uncorrelated. 


Table 17: Correlated and uncorrelated systematic uncertainties. The / symbol indicates that 
the uncertainty has been taken into account in the analysis, but it is not correlated with any of 
the other analyses. The EZl symbol indicates that the uncertainty has been taken into account 
and that it is correlated with the other analysis that have a EZl sign as well. A missing symbol 
indicates that this uncertainty is not relevant for this analysis charmel. 



Single 

Inclusive 

Multiple 

All-had. 

All-had. 

H —> 77 
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7 Results 


7 Results 


No significant deviation from the SM prediction is observed. The expected limits of the indi¬ 
vidual analysis channels at a 95% confidence level (CL) are displayed in Fig. for exclusive 
decays of the T quark to tH, tZ, and bW. This figure also shows the result of the combination, 
where only the non-overlapping part of the individual analyses are combined, as discussed in 
Section m 





Figure 7: Expected limits at 95% CL of the individual analyses in comparison to the combina¬ 
tion for exclusive decays of the T quark to tH, tZ, and bW. 


The observed limits and the expected one and two standard deviation uncertainties are dis¬ 
played in Fig. I^for exclusive T quark decays. 


The lower limits on the mass of the T quark are obtained by determining the intersection be¬ 
tween expected (observed) limits with the theoretical prediction, based on the cross section 
versus T quark mass distributions shown in Fig.[^ The results are visualized graphically in the 
triangular plane of branching fractions in Fig. 1^ The numerical upper limits on the T quark 
production cross section are given in Table 18 for a full range of branching fractions and the 
numerical results of the limits on the mass of the T quark are given in Table 19 A different 


visualization of the mass limits is presented in Fig. 10 


Depending on the assumed branching fractions, the expected limits lie between 790 and 890 GeV, 
while the observed limits are in a range between 720 and 920 GeV. In much of the triangular 
plane of branching fractions these are the most stringent limits on T quark pair production to 
date. 
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Figure 8: Observed and expected Bayesian upper limits at 95% CL on the T quark production 
cross section for exclusive T quark decays to tH, tZ, and bW. The green (inner) and yellow 
(outer) bands show the 1(7 (2cr) uncertainty ranges in the expected limits, respectively The 
dashed line shows the prediction of the theory 
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7 Results 


Table 18: Branching fractions (first three columns) and the observed and expected upper limits 
on the T quark production cross section at 95% CL for different values of the T quark mass. 
The expected limits are quoted with their corresponding uncertainties, for different branching 
fractions hypotheses. The cross section limits are given in units of pb. 
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Table 19: Lower limits on the mass of the T quark at 95% CL, for different combinations of T 
quark branching fractions. The Itr uncertainty range on the expected limits are given as well. 


B{tH) 

BitZ) 

e(bw) 
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8 Summary 

A search for pair production of vector-like T quarks of charge 2/3 has been performed. In most 
models the hypothetical T quark has three decay modes: T —> tH, T —> tZ, and T —> bW. The 
following five distinct topologies have been investigated: inclusive lepton events covering all 
possible decay modes, single-lepton events optimized to find T —> bW decays, all-hadronic 
events optimized either for T —> tH or T —> bW decays, and events containing a Higgs boson 
decaying to a pair of photons. 

Data and SM background expectations are found to be in agreement. Upper limits on the 
production cross sections of vector-like T quarks are set. The expected 95% CL lower mass 
limits are between 790 and 890 GeV depending on the branching fraction of the T quark. For 
a branching fraction of 13 (tH) = 100% an expected (observed) limit of 840 (770) GeV is found. 
For 13(tZ) = 100% the expected (observed) limit is 830 (790) GeV and for 13(bW) = 100% the 
limit is 890 (920) GeV. These are among the strongest limits on vector-like T quarks obtained to 
date. 
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